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Abstract 

It is shown on the basis of scahng arguments that a disordered 
interface between two elastic soHds will quite generally exhibit static 
and "dry friction" (i.e., kinetic friction which does not vanish as the 
sliding velocity approaches zero), because of Tomlinson model insta- 
bilities that occur for small length scale asperities. This provides a 
possible explanation for why static and "dry" friction are virtually 
always observed, and superlubricity almost never occurs. 

1 Introduction 

A few years ago, Muser and Robbinsjl] proposed that static friction be- 
tween clean crystalline solid surfaces, whose crystal axes were rotated by an 
arbitrary angle with respect to each other so that the surfaces are incom- 
mensurate, would be zero. This is commonly known as " superlubricity" [2] . 
In order to explain the fact that superlubricity is rarely observed, they then 
proposed that the likely source of static friction for such surfaces is the exis- 
tence of mobile dirt molecules, always present at any interface. The mobile 
molecules seek out local energy minima at the interface, which results in the 
surfaces being pinned with respect to each other, i.e., there is static friction. 



Actually, there can also be superlubricity for smooth disordered surfaces. 
In Refs. [21 S] it was argued that a fiat disordered interface between two 
macroscopic size surfaces which do not interact chemically will exhibit effec- 
tively no static friction for interface interaction per unit interface area small 
compared to the shear elastic constant and high friction for interface interac- 
tion above this value. This frictionless regime is known as the weak pinning 
regime and the large friction regime is known as the strong pinning regime. 
In Ref. [3], I argued that the existence of micron length scale asperities at 
an interface can result in both slow speed kinetic (often referred to as "dry 
friction" ) and static friction, resulting from the fact that these asperities can 
occur in multistable equilibrium configurations, which result in " dry friction" 
as a consequence of the Tomlinson model [5] and static friction as well. The 
requirement for asperities being multistable is that the force on an asperity 
on one surface due to its interaction with a second surface dominate over the 
elastic forces resulting from its distortion, as a consequence of its interaction 
with the second surface. Since the interaction of an asperity with the second 
surface varies on atomic length scales, the distance over which the asperity 
distorts is negligibly small compared with its size (about a micron). Since 
micron length scales are still large compared with atomic spacings, however, 
the resulting static and dry friction is still likely to be quite small, indicating 
that the Muser-Robbins picture is still basically correct. 

It is argued here that sub-micron scale roughness can have a crucial effect 
on static friction between two three dimensional elastic solids that interact at 
their interface, primarily through inter-atomic hard core interaction, which 
will be the case for sufficiently high load [6J. The basic idea is that the 
multiscale nature of the surface roughness results in the load being supported 
by a very small percentage of the atoms that would be in contact with the 
second surface, if the surface were perfectly fiat. Consequently, the forces 
acting on these atoms by the second surface will dominate over the elastic 
restoring forces. Therefore, the shorter length scale asperities will be quite 
likely to exhibit multistability. This mechanism was previously proposed 
as a possible mechanism for explaining the ability of hard surface coatings 
to function as good lubricants [7]. Plasticity (which was not discussed in 
Ref. [7J but is discussed here) will likely play an important role for small 
length scale asperities as well. It will be argued that plastic deformation 
under the high degree of stress acting on the small length scale asperities 
will enhance the tendency for multistability to occur for these asperities, 
resulting in an increased tendency for the occurrence of static and "dry" 
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friction. Recent work on contact mechanics [HI E] niake it possible to provide 
a more complete treatment of the scaling theory of Ref . [7] , allowing a more 
numerical study of the conditions under which static friction can arise as 
a consequence of multiscale roughness, without the need to postulate that 
static and dry friction can only occur when there are mobile ion dirt particles 
present at the interfacepQ, at least for surfaces that can be modeled as self- 
affine over a few length scales. 

The next two sections summarize the results of Ref. [7j and make some 
improvements and corrections. In section IV, the scaling theory of Ref. [7] 
is expressed in the language of Persson's theory of contact mechanics[9], and 
the results of this theory are used to obtain information on the conditions 
under which large static and "dry" friction are likely to occur. 

2 Static Friction due to Hard Core interac- 
tions 

As in Ref. [7], we consider two surfaces in contact which are disordered so 
that those atoms which are in contact are randomly distributed over the 
interface. In order to simplify the discussion, we will consider as our model 
for the interface, a rough elastic surface in contact with a perfectly atomically 
flat rigid surface. This has been argued in the past to be equivalent to the 
more correct model of two rough elastic surfaces in contact, for the contact 
mechanics problem, in which there are only forces normal to the surface [T0| 
m [8]. While it is not a rigorously correct representation of the problem 
when there are frictional (i.e., shear) stresses present as well, it should give 
correct orders of magnitude for the scaling treatment considered on Ref. [7] 
and here. Let us also assume that the atoms in contact at the interface 
interact only with hard core interactions. This could occur either because 
the surface atoms are chemically inert and there is negligible adhesion, or 
because the surfaces are being pushed together with a sufficient load so that 
the hard core interactions dominate. Let a denote the load or normal force 
per unit interface the mean atomic spacing and c, the fraction of 

the surface atoms of one surface that are in contact with the second surface. 
Then, each of the atoms in contact must contribute on the average to the 
normal force a force of order aa? /c. Since the force due to the hard core 
interaction between a pair of atoms acts along the line joining the atoms, for 
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Figure 1: This figure illustrates how the hard core interaction between a 
pair of atoms, one belonging to each of the surfaces in contact, can both 
support the load and give rise to static friction between the surfaces. Since 
the force F between the pair of atoms can have both a component normal 
to the interface, F^, which contributes to the normal force supporting the 
load, and a component along the interface F^,, the mean value of must be 
proportional to the mean value of F^. 

most relative positions of the atoms, it has a component along the interface, 
as illustrated in Fig. 1. 

In the strong pinning regime, each surface atom will sink into an interface 
potential minimum at the expense of the elastic forces holding it in place. 
Such a minimum will generally occur at an interstitial region on the second 
surface. If we now attempt to slide the surfaces relative to each other, each 
of the atoms in contact with the second surface will now, as it gets pulled 
out of its potential minimum, exert a component of its hard core interaction 
with the second surface parallel to the interface directed so as to oppose the 
attempted sliding motion. This is identified with the static friction. Since 
each of these atoms must also provide a component aa^/c normal to the 
interface on the average, it is clear the static friction is proportional to the 
load. This is the case because the contribution to the load and the static 
friction for each pair of atom in contact is provided by the same hard core 
force acting between the atoms. The proportionality constant Hs is identified 
with the coefficient of static friction which is not too much smaller than 1. 
This accounts for Amonton's law without the need to assume that the friction 
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is proportional to an ill-defined area of real contact. In the weak pinning 
limit, the component along the interface of the hard core force is random, 
and hence, for an infinite interface area and hence an infinite number of 
interface atoms, the components along the interface of the hard core forces 
cancel, resulting in effectively no static friction in the thermodynamic or 
macroscopic solid limit. 

3 Effects of Roughness (i.e., Asperities on 
Several Length Scales) 

For two perfectly flat surfaces in contact pushed together with sufficiently 
weak load, each containing N surface atoms, the mean force of static friction 
will be the sum of the components of these hard core forces parallel to the 
surface, which have random magnitudes and directions. The net static fric- 
tion, which is the sum of these components, is proportional to N^^"^ because 
they do not add together coherently. This is known as the weak pinning 
regime [m [121 [13]. It was shown in Ref. [IH [121 [131 [3] that for sufficiently 
weak interaction across the interface (which will be true at low loads), the 
system will be in this regime. (Although mica is probably the only mate- 
rial with complete atomic scale flatness, bare mica surfaces actually exhibit 
relatively high friction because they interact with each other quite strongly.) 
At higher loads, the atoms from each surface can be pushed sufficiently far 
into the regions between atoms on the second surface (resulting in the atoms 
from the second surface being pushed apart), as this will minimize the re- 
pulsive energy. When we attempt to slide the surfaces relative to each other, 
there will be a net component of hard core repulsive forces along the surface 
opposing the sliding, i.e., static friction. In this case, the net force of static 
friction will be proportional to N. For models for surfaces with single length 
scale roughness, like the Greenwood- Williamson model [H], since there are 
of the order of 10* atoms at an interface between two micron size asperities, 
since N^/"^ is 10^, the static friction is a factor of 10^ smaller in the weak 
pinning than in the strong pinning regime. From Ref. we find that for 
flat surfaces, the load per unit area a below which the interface is in the 
weak pinning regime is given by a ~ c^^'^K^ where c is the fraction of surface 
atoms in contact with the substrate and K is the shear elastic constant. Since 
K ^ l{)^^N/m? for most solid materials, unless c is extremely small, we find 
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on the basis of this argument that any flat surfaces in contact that interact 
with only hard core interactions will be in the weak pinning regime. 

Of course, no surfaces are perfectly smooth, and in fact, there can be 
roughness on several length scales. Let us assume then that there are rim 
orders of length scales, which we represent as follows: We divide the surface, 
which is defined to be at the length scale, somewhat arbitrarily into 
M„^_i regions, which we define to be the asperities at the — 1** order 
length scale, a fraction c„„_i of which are in contact with the second surface. 
An asperity at this length scale is considered to be in contact with the second 
surface if some part of it is in contact with the second surface. Let N represent 
the number of atomic size regions (i.e., 2-dimensional unit cells) contained 
in the projection of the rough surface onto a perfectly fiat surface. In the 
discussion in this paragraph, the asperities are treated as rigid " hills" on the 
surface. Effects of distortion of the asperities are treated later in this section. 
As mentioned above, at each pair of asperities in "contact" the contact is 
likely to only occur at selected isolated regions, which we may refer to as 
rim — l"** order asperities. In fact, we will divide the surface of each — 1** 
order asperity into M„^_2 rim ~ 2"^^ order regions, somewhat arbitrarily (as 
we did for the rim — 1** order asperities), a fraction c„„_2 are in contact, 
again, contact meaning that there are high regions within this region which 
are in contact. The interface between a pair of rim — 2 order asperities can be 
divided up into rim — 3 order asperities as well, M„^_3 asperities, a fraction 
Cnm-3 of which are in contact, etc., until we have done rim sub-levels of this 
sub-division. The surfaces can never be truly self-affinej9], however, because 
when we reach atomic dimensions at n=0 order, this sub-division into smaller 
and smaller length scales terminates. The area of contact of a 0*'' level (i.e., 

smallest) asperity will contain of the order of A^n = tt — rf ti 

atoms, a fraction Ca of which are in contact. It is these atoms at the zeroth 
(i.e., the final) order of asperities which support the load. 

If the n*^ order set of asperities are in the weak pinning regime, the static 
friction acting on it is reduced by a factor (c„M„)^/^, because by the above 
arguments, the static friction forces from these asperities act incoherently, 
and if the atoms at the interfaces of O*'^ order mini-asperities in contact are in 
the weak pinning regime, the static friction is reduced from the strong pinning 
regime value (i.e., fig not too much smaller than 1) by a factor (caA'o)^/^. This 
the mechanism was proposed as a possible way to explain why coatings of 
stiff materials are good lubricants [T5]. 
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The arguments in the above paragraphs are only correct if each asperity 
is completely rigid, as we have assumed that the elastic forces which oppose 
displacements of the points of contact with the substrate resulting from the 
forces that the substrate exerts on them are due to the bulk solid's elasticity. 
Since each asperity has some height, however, it can distort so as to move its 
interface with the substrate closer to its potential minimum, even without 
distorting the bulk solid. Thus, whereas a particular substrate with com- 
pletely rigid asperities might be in the weak pinning regime, if the asperities 
are able to distort by a sufficient amount, it might be in the strong pinning 
regime. This may explain why it is that although the estimates given ear- 
her in this section indicate that most solids should be in the weak pinning 
regime, this is not consistent with the magnitudes of the friction coefficients 
that are observed for most solids. 

Let us now consider the distortions of the asperities that occur in response 
to the substrate potential. In our model, there exists a bunch of smallest 
asperities (which will be considered the lowest or zeroth level) which are 
in contact with the substrate. There are several groups of these that are 
assumed to be attached to a bunch of larger next or first order asperities. 
Groups of these first order asperities are then attached to larger asperities 
which are known as second order asperities. This hierarchy continues until 
we reach an "asperity" of width equal to that of the whole interface. This 
hierarchy of asperities is illustrated schematically in Fig. 2. Note that this 
figure is only a schematic representation of the problem; it is not being 
suggested here that the asperities have a rectalinear shape. Possible effects 
of asperity shape are discussed in appendix A. 

Consider the zeroth, the lowest order (i.e., the smallest), asperity. Let 
it have a height of order Lq and a width of order Lq. To find its distortion 
resulting from the sum of the substrate potential energies of all of the atoms 
of the asperity which are in contact with the substrate, we must minimize the 
sum of its elastic and substrate potential energies. The substrate potential 
energy is given by Vo{c]/'^Lo/a)fo{A.xo/a), where Vq is the amplitude of the 
interaction of a single atom with the substrate, resulting primarily from hard 
core repulsions between the atoms, Axq is the amount that the surface in 
contact with the substrate slides under the infiuence of the substrate potential 
as the asperity distorts while all higher level asperities remain in an arbitrary 
rigid configuration and /o(Axo/a) is a function of order unity which gives the 
variation of the substrate potential with Axq for fixed, undistorted asperities 
of higher order (i.e., larger size in the present context). (Clearly, each of the 
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Asperity Hierarchy 



Figure 2: This is a schematic illustration of the asperity hierarchy on the 
top surface sliding on a flat substrate (i.e., the bottom block). (Real asper- 
ities have arbitrary shapes, as opposed to the square shapes shown in this 
schematic representation.) Each asperity of a given order has a number of 
(smaller) asperities of one order lower on its surface. In turn, each of these 
lower order asperities has a number of (smaller) asperities of one order lower. 
This continues until we reach the zeroth order asperity, whose surface consists 
of atoms, although only three orders of asperities are illustrated here. 

zeroth order asperities has a different function; /o denotes a generic function 
describing the interface potential energy for a typical zeroth order asperity.) 
Each function clearly must possess multiple minima. We are assuming here 
that the surface of the asperity in contact with the substrate is in the weak 
pinning limit. The factor (cy^Lo/a), which is of the order of the square root 
of the number of atoms in this surface expresses this fact. If the surface 
of the asperity in contact with the substrate is in the strong pinning limit 
instead, this factor will be replaced by Ca{Lo/a)^, the number of atoms at 
the interface. [The factor of c^^ and Ca, respectively in these expressions 
were inadvertently not included in Ref. [7].] Treating this asperity as an 
elastic three dimensional solid in contact with the substrate we find from 
the discussion in Ref. [4j that the interface between the substrate and this 
asperity is in the weak pinning limit if o"o < cl^'^K, where Uo is the mean 
load per unit interface area supported by this asperity and is the fraction 
of the surface atoms of this asperity that are in contact with the substrate. 
Assume that a fraction cq of the zeroth order asperities have atoms belonging 
to them in contact with the substrate. Let ci represent the fraction of next 
order (first order) asperities whose zeroth order asperities are in contact with 
the substrate, C2, the fraction of second order asperities whose first order 
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asperities have their zeroth order asperities in contact with the substrate, 
etc., up to order. Then ctq = cr/(coCiC2...c„^„i), where a is the load 
per unit apparent area of the surface of the whole solid. Then, we conclude 
that the criterion for the atoms at the interface between the zeroth order 
asperity and the substrate to be in the weak pinning regime is that a < 
{c]/'^CQCiC2-..Cn^-i)K . We see from this inequality that the more fractal the 
surface is, the more difficult it is for the zeroth order asperity to be in the 
weak pinning regime. The cost in elastic energy due to the shear distortion 
of the asperity can be determined by the following scaling argument: The 
elastic energy density for shear distortion of the asperity is proportional to 
[duxjdzf' ■, where represents the local displacement due to the distortion, 
the x-direction is along the interface and the z-direction is normal to it. 
The Ux must scale with Axq and the dependence of Ux on z has a length 
scale Lq. Thus the elastic strain energy of the asperity is of the order of 
(l/2)LoI/Qii'(Axo/I/o)^, where K is the shear elastic constant and (Axo/I/q) 
is the average shear strain and Lq is the mean width of the asperity. Using 
a perhaps more correct non-rectangular shape for the asperities is shown in 
appendix A to only produce a correction of order unity. Minimizing the sum 
of these expressions for elastic and substrate potential energy, we obtain 

Axo/a ^ {y^lKa''yl\L'jU)r^{^x^la) (1) 

since /q, the derivative of /o with respect to its argument, it is of order one, 
from the definition of /q. Let us follow a line of reasoning like that of Ref. 
[5], a modified version of which is given in Ref. [7]. For iy^j Kc^^c^J^'^iL'^j L^) 
below a certain value of order one, for small Vo/i^a^, Eq. (1) can have only 
one solution for Axq. The reason for this is illustrated in Fig. 3. Under 
such circumstances the average kinetic friction, in the limit as the sliding 
velocity approaches zero, is zero. For a surface with an infinite number of 
asperities, distributed uniformly in space, it was shown in Ref. [5j that the 
static friction is zero as well. A modified version of this argument, which 
points out that for a surface with a finite number of asperities the static 
friction is non-zero, but smaller by a factor of a/Lg compared to what it would 
be if the contributions of the asperities to static friction acted coherently is 
provided in appendix B of Ref. [7]. If this asperity is in the strong pinning 
limit instead, we replace the factor of cy^(Lo/a) by Ca{LQ/aY to account for 
this and as a result, the factor c]J'^{L'q/ Lq) gets replaced by CaiLya), which 
could make the asperity satisfy the criterion for multistability more easily, 
and consequently, the friction from these asperities will no longer be reduced 
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by the factor a/L^. For a load per unit area a, assumed to be primarily 
due to hard core interactions, we may assume Vq ~ aa^/c, where c is the 
fraction of the surface atoms which are in contact with the substrate. By 
the above arguments, c = CaCoCiC2...Cn,„-i. Then, we see that the criterion 
for the zeroth order asperity to be multistable is a > CQCiC2...Cn^-iK . If the 
criterion for weak pinning for the zeroth order asperity surface is not satisfied, 
the criterion for monostability of this asperity gets changed from the above 
inequality to {Vq/ Ko?)ca{L'Q/ a) < 1, which is more difficult to satisfy since 
Lq/o can be considerably greater than 1. 

At the next level, we have an asperity surface in contact with the sub- 
strate which consists of a collection of the lowest level (i.e., the smallest) 
asperities discussed in the previous paragraph. Assuming this asperity to 
be in the weak pinning regime, the potential of interaction with the sub- 
strate, which is the sum of all of the interactions of the substrate with 
the lowest order asperities, which cover a first order asperity, is of order 
Vo(caCo)^/^(Lo/a)(Li/Lo)/i(Aa;i/a). Here Li represents the width of this 
order asperity, Axi represents a displacement of the lower surface of this 
level asperity for fixed (i.e., un-distorted) configurations of all higher order 
asperities, and /i denotes one of the functions which describes the interface 
potential energy of one of the first order asperities. It has at least one min- 
imum and runs over a range of magnitude one as its argument runs over a 
range of order one. [Again, the factor of (caCo)^''^ and c^Co, respectively in 
these expressions were inadvertently not included in Ref. |7].]The elastic en- 
ergy is of the order of (l/2)L'^Lfi^'(Axi/L'^)^, by the argument given above 
Eq. (1), where L'-^ is the height of the body of the first order asperity, which 
is assumed to be much larger than Lq. Minimizing the sum of these two 
energies, we obtain 

Axi/a ^ {Vo/Ka^){caC^Y'\L'jL,)f[{^x,/a). (2) 

Again, we conclude, based on the arguments presented in Ref. [5], that 
the static friction is reduced by a factor of Lq/Li below what it would be 
if the contributions to the static friction from each of the mini-asperities at 
this level acted coherently. If the zeroth order asperities attached to this first 
order asperity are in the strong pinning regime, the factor of {caCQY^'^Li/ Lq in 
the equation for the interaction of this asperity with the substrate is replaced 
by CaCo(Li/Lo)^, and hence, the right hand side of Eq. (2) has the factor 
{caCQY^'^L'^l Li replaced by CaCQL[/LQ, which can make the solutions to this 
equation for Axi multistable. 
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Figure 3: This figure illustrates the solution of Eqs. (1), (2) and (3), for 
Axq, Axi and Ax„, respectively. f'{x) is a schematic illustration of the 
functions /q, f[ and and x denotes Axq, Axi or Axn, respectively. Lines 
A and B represent the line y = {Ka^ /Vo){Ln/ L'^)x, for Ka^ /Vo){Ln/ L'^) < 1 
and Ka^ /Vo){Ln/ L'^) > 1, respectively. For the situation illustrated by line 
A, there are multiple solutions (i.e., multi-stability), while for the situation 
illustrated by line B, there is only one (mono-stability). 

Continuing this procedure, we find that the displacement of the n^^ level 
mini-asperity is found by solving 

Axja ^ {Vo/Ka^){ca...Cn-iy^\L'jL^)f:^{AxJa), (3) 

where L„ and L'^ are the width and height of the body of the n*'* level as- 
perity. [The factor of (ca...c„_i)^/^, respectively in these expressions were 
inadvertently not included in Ref. [7]. The general conclusions of Ref. 
|7j, however, are not changed by this omission.] If (Vo/-ft'a^)(ca...c„_i)^/^ = 
a/[K{cnCn+i---Cn^y^'^] > 1, and L'^/Ln ~ 1 for all n, asperities of all orders 
will be multistable, implying the occurrence of large static friction. If the 
condition given earlier for strong pinning at the zeroth order asperity inter- 
face, namely a > cl^^CoCi...Cn^K = {c/cl^'^)K is satisfied, the condition for 
multistability on all levels, a > {cnCn+i-.-Cn^Y^'^K is certainly satisfied for 
the zeroth order asperities, but may break down at some order n. Since more 
of the c's will appear in the product of the c's in this expression for lower or- 
der than higher order asperities, we conclude that the lower order asperities 
are more likely to be multistable than the higher order ones. If the n^^ order 
asperity is in the strong pinning regime, the factor (ca...Cn_i)^/^(L^/Ln) in 
Eq. (3) gets replaced by (ca...c„_i)(L^/a), which leads to the condition for 
multistability: a > (c„c„+i...c„„)(L^/a)iir. 

At least in the limit in which the stress acting on a given length scale 
asperity resulting from the load that it carries is small compared to Young's 
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modulus, recent theories of contact mechanics for surfaces with multi-length- 
scale roughness [H [H] show that the area of contact on all length scales (even 
the smallest) is approximately proportional to the applied load. This implies 
that although the smallest length scale asperities will certainly flatten out as 
the load increases (which means in the language of the present manuscript 
that a larger fraction of the subasperities on the surface of this asperity will 
come in contact with the second surface, and a similar picture holds for each 
of its subasperities). Consequently, on the basis of these theories, the amount 
of load carried by each of the atoms which are in contact with the second 
surface will not decrease as a consequence of more of these atoms being in 
contact with the second surface. 

4 Application of Persson's Contact Mechan- 
ics Theory 

Persson has developed an analytic contact mechanics theory for self-affine 
surfaces[9]. In this section, a connection will be made between the parame- 
ters used here and the contact area on a given length scale A, A{() used in 
Ref. [9], were ( = L/X, where L is the largest length scale of the surface. On 
a length scale A = there will be (L„^/L„^_i)^ asperities of lateral size 
Lum-i, a fraction c„^_i of which are in contact with the substrate. It is clear 
that Ln„-, denotes the largest length scale, which represents the dimensions 
of the actual surface. On this surface, there will be c„„_2(-C/n™-i/-^n„-2)^ 
contacting asperities of size L„^_2, and on each of these asperities, there will 
be Cn^-3{Ln^_2/ Ln^_3y coutactiug asperities of size Iv„^_3. As we continue, 
we reach the Ln+i scale asperities, each of which contains Cn{Ln+i/ Ln)^ con- 
tacting asperities of linear dimension L„. If we were to stop at this length 
scale and assume that the area of contact of each of these L„ length scale 
asperities is smooth (because we are imagining that our "measuring instru- 
ments" cannot see such small length scales), and the total area of contact of 
the surface with the substrate will be given by the product of the above num- 
bers of contacting asperities at each length scale and (the cross-sectional 
area of an n^^ scale asperity), or 

P{Cn) = A(Cn)M(C„,J = 
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^2 

Cn^-?.{Lnm-2/ Ln^-i) ...Cn{Ln+l/ Ln) (-^) = 

Cn„-lC„„_2C„^-3---Cn+lCn, (4) 

where we have used the fact that y4(L„^) = Iv^^^ , and hence, the quantity 
P(Cn) = A{.Ln) / A{Ln^) is a quantity used in Persson's theory p], where the 
magnification = The largest length scale -Ln„, which will often 

be denoted by L, is the length of sliding solid, and A = L"^ is its nominal 
area. Using this relationship between Persson's -P(Cn) and the quantity c„ 
used in Ref. [7] and here, we can write the criterion for multistability of the 
^th Qj^^gj^ asperity discussed under Eq. (3) of the last section as 

^ > {a/L'^)P{Cn)K = (5) 



assuming strong pinning at all asperity interfaces and 

LL A 



a > {LjL'^)[P{Cn)P{QY''K = ^l^M^M^K, (6) 



n 



where (a = Ln^/a, assuming weak pinning. From Eqs. (5) and (6), it is clear 
that for sufficiently small ^(Cn), the asperities will tend to be multistable, 
resulting in nonzero static and dry friction. By writing the criteria for multi- 
stabihty of the asperities in terms of Persson's -P(Cn), we can use the results 
of his theory of contact mechanics to determine the load dependence of these 
criteria. 

Let us now use Persson's theory to determine the contact area per asperity 
at a particular length scale (call it the tt}^ length scale). Since the number 
of rim — 1 level asperities in contact is c„^_i(L„^/L„,^_i)^, and the number 
of rim — 2 level asperities residing on one of the rim — 1 level asperities in 
contact is c„„-2(-^n„-i/-^n,„-2)^5 etc., we conclude that the number of n*'^ 
level asperities in contact is 

Ln^~zf ■■■Cn{Ln+l/ Lnf , (7) 

which is equal to 

Cn^-lCn^-2Cnm-^--Cn{Ln^/Lnf = A(L„)/L^^. (8) 
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Using the fact that the contact area per n order asperity is equal to the 
nominal area of the interface A{Ln^) divided by the above expression for 
the number of n^^ level asperities in contact, we obtain for the mean contact 
area per asperity at the n^^ level L^, independent of the load. If we stop at 
^th Qj^^gp jj-^ Persson's treatment [9] , we take the surfaces at smaller length 
scales to be perfectly flat. It is interesting to note that if we calculate the 
area per n^^ order asperity using Archard's theory of contact mechanics for 
fractal surfaces p^. we find that it is proportional to the load to a power 
which decreases as n decreases, implying that even if the surface is assumed 
to be perfectly elastic, the smaller length scale asperities will fiatten out 
at sufficiently high load, provided that we put a small length scale cut-off 
in his theory. In reality, there must also be a small enough length scale in 
Persson's theory as well in which asperities at all smaller length scales will get 
completely flattened out, if the solid is treated as being completely elastic, 
but of course long before that point, we must deal with plasticity. 

Persson shows using his theory of contact mechanics [Q] the not unexpected 
result that as the magniflcation parameter (n = L^^/Ln approaches inflnity 
^eKCn)M5 where Aei{Cn) denotes the area of contact associated with elastic 
deformation of asperities, approaches zero. Then EyApi^Cn) = ctqAq, where 
Api{Cn) is the contact area associated with plastic deformation of asperities 
under load and Ey denotes the hardness of the material (i.e., the compres- 
sional stress at which plastic deformation takes place). Some of Persson's ar- 
guments are summarized and some steps are fllled in appendix B. Generally, 
when lowest length scale asperities have failed plastically under compression, 
if one attempts to shear the asperity by sliding the surfaces relative to each 
other, it will also have failed plastically in shear motion|18j, and hence, ei- 
ther the resulting shear stress due to shearing of the asperity will be nearly 
independent of the shear strain or it will increase with shear strain , but 
much more slowly, as a function of shear strain, than it would if its response 
were elastic. In the latter case, the tendency towards multist ability will be 
greatly increased over what it would be if the shear response were elastic, 
implying an increase in the dry and static friction. In the former case, the 
asperity would shear without any increased restoring force, most likely indi- 
cating that it will break as it shears relative to the substrate. This indicates 
that the system has switched over to a regime in which a good part of the 
kinetic friction is due to wear (i.e., the breaking of the smallest length scale 
asperities). The static friction will likely go to zero if this were a completely 
correct model for plastic flow, because in such there is no energy cost 
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to shear the asperity. 

As mentioned earher, Persson has demonstrated that as (n approaches 
infinity, all parts of the contact area exhibit plasticity. (The argument is 
discussed in appendix B of this article). It is important, however, to esti- 
mate the value of („, for which plastic deformation of asperities on the ra*'* 
scale becomes important. The rough criterion that we will use to estimate 
this is when the mean value of the compressional stress over an n'^'^ order 
asperity becomes comparable to Ey- One may determine when plasticity at 
the smaller length scale asperities sets in using Persson's theory, assuming 
that the distortions on all length scales are elastic until this point. Persson 
finds in the elastic regime the following expression for the contact area at n*'' 
scale magnification: 

^ . P(c„) . l::(i_i^(^£)v.e-... (9) 

A qoiioE 7iH 

where H = 3 — Df, where -D/ is the fractal dimension of the surface, go = 
271 /L and ho is the rms surface height fiuctuation, and u is Poisson's ratio. 
Since we must have 

<a> |^„A(C„) = aA, (10) 
< > Icn = '^/P{Cn), and hence, from Eq. (9), we have 

< ^ > kn _ gp/^o . TiH .i/2>(i-j/) 



E 4(1 -z/2)'l -iJ' 



C""^- (11) 



When < cr > becomes comparable to Ey, one must consider plasticity. 
Persson considers H=0.8, a value typical for pavement and go^o equal to 0.001 
and 0.01. [Smaller values of H represent rougher surfaces, because by the 
definition of a self-affine surface, reducing the length scale along the surface 
by a factor A (with A < 1) reduces the height fiuctuations only by a factor A^, 
and since H < 1, larger height fiuctuations occur over smaller lateral length 
scales.] Using these numbers in the above expression for < cr > |^^, we find 
that for L/Ln = 10^, typical for atomic level asperities, < a > \(„/E ^ 0.04 
for qoho = 0.001 and 0.4 for go/io = 0.01. Ey/E=0.1 for glass, 0.0013 for 
steel and 0.0125 for hard plastic. Thus, atomic scale asperities are clearly 
in the plastic regime. For diamond films, Ey, the hardness, is about 100 
GPa[19j, and since E is 1080GPa[20j, Ey/E ^ 0.09. Other diamond-like 
carbon films with large sp^ (i.e., tetrahedral a) bond content also exhibit 
Ey/E of comparable magnitude. We see here the important role played by 
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the relatively large values of Ey/E that exist for these materials. On the 
basis of the above arguments, we conclude that it is possible for sufficiently 
smooth surfaces that for diamond-like carbon films, asperities on all length 
scales that exist for a crystalline material might be elastic. 

At this point, let us re-examine the argument presented in Ref. |7j, that 
the large values of Young's modulus of stiff coatings, such as diamond-like 
carbon films, can still account for their excellent lubricating properties, mak- 
ing use of Eqs. (5), (6) and (9). Combining Eqs. (5) and (6) with Eq. (9), 
we obtain for the condition for the n*^ order asperities to be multistable 

E^ a 4(1 -z/2) 4-/7^ ' ^ 

if the interfaces between each level asperity less than or equal to n*'^ order 
are in the strong pinning regime (defined in the previous two sections), or 

E^ L^A{l-u^)^l-H^ ^^"^'^^ ' ^^^^^ 

if the interfaces between each level asperity less than or equal to n*'* order 
are in the weak pinning regime. Since the n=0 order asperities are the ones 
most likely to be multistable, let us apply Eqs. (12a) and (12b) to the n=0 
case, for which they reduce to 

K L'q qphp TvH .i/2.(i-H) ..ON 
E ^ a 4{l-u^)h-H^ ' ^^'^"'^ 

respectively. The condition for weak pinning at an n=0 asperity interface 
found in the last section is 

a < c'J'{cpC,...Cr.^.,)K, (14) 

which when written in terms of Persson's notation [H] is 

a < c-'/^^Mk. (15) 
A 



and 
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Substituting for A{(a)/A using Eq. (9), this condition can be written as 



In order for the zeroth order asperities to be muhistable, and at the same 
time have their interfaces with the substrate in the weak pinning hmit, we 
must have c^^ < L'q/Lq, which is not difficult to satisfy unless the zeroth 
order asperities are all extremely short and fat and/or have a large fraction 
of the atoms at their interfaces in contact with the second surface. 

5 Conclusions 

It has been shown, using a scaling theory of friction for surfaces with mul- 
tiscale roughness^ combined with Persson's theory of contact mechanics, 
that, in contrast to surfaces with only single scale roughness which appear to 
always exhibit superlubricity, unless one postulates the existence of mobile 
dirt particles at the interface, surfaces with multiscale roughness will al- 
most always exhibit static and dry friction, characteristic of almost all solid 
surfaces. The static and dry friction come about because asperities at the 
smallest length scales will be multistable, because they support all of the 
load if they remain elastic, or will exhibit plastic failure. It is argued that 
plastic failure will increase the tendency towards the asperities being multi- 
stable. Numerical finite element calculations used recently to study contact 
mechanics [S] can also be applied to this problem. Since they appear to yield 
results for the contact mechanics problem which are qualitatively similar to 
those of Ref. [9], however, the results for the multistability (i.e. Tomlinson 
model[5]) problem should not differ significantly from those presented here. 
There is one small difference, however, which is worth considering. Whereas 
Ref. [9] finds that A{(n)/A is a linear function of a in the elastic regime, Ref. 
[H| reports that this function is a slightly sub-linear function of a. As already 
suggested in Ref. |7J , this would imply that as a increases, the interface could 
possibly switch from having all asperities in the monostable regime, resulting 
in low friction, to the multistable regime, resulting in high friction, when a 
becomes sufficiently large. 




(16) 



17 



6 Appendix A: Effects of Asperity Shape on 
Friction 



In section III, we assumed that the actual shape of an asperity would have 
an insignificant effect on the multiscale Tomlinson model that has been de- 
veloped here. In this section, we will examine whether this is likely to be a 
reasonable assumption. Asperities on all length scales are likely to initially 
be peaked, but as they are pressed together they fiatten out. The amount 
of flattening of a given length scale asperity depends on the amount of load 
that asperity carries, which itself depends on the degree of compression of 
the asperity, as is known to happen for single length scale asperities pLQj- The 
amount of load carried by an asperity depends on the number of atoms in 
contact with the substrate, as we have seen earlier that each atom carries a 
load of the order of Pa^/c. Typically, near the peak of an asperity at any 
length scale, the proflle can be represented to a good approximation by a 
parabaloid whose cross-sectional area at a distance z from apex is propor- 
tional to z[TO]. It is not a bad approximation to assume that we can use this 
result for an asperity which gets compressed because it is in contact with the 
substrate. Then let us assume that a compressed asperity has a flat peak at 
a minimum value of z (-Zmm), whose area is proportional to Zmin- Then, using 
the parabaloid approximation for the asperity, we have vrr^ 2nRz, let us 
imagine slicing each such asperity into slices of thickness Az, each of whose 
area is proportional to z. Then, the shear stress between slices n and n+1 is 
equal to 

a a 
following the method for taking the continuum limit of the harmonic ap- 
proximation treatment of a discrete lattice discussed in Ref . [1^ . This must 
be zero when the asperity is in equilibrium. The continuum limit of this 
equation is 

d d{zu) 



dz dz 

whose solution is 



0, (2A) 



az z 

where c is a constant. Since aX z = Zmin, we require that K{'Kr^){du/ dz)\z=z„,i„ 
crvrr^, where is the contact radius, and a is the shear stress acting on the 
contact area of the asperity, because this is the condition for equilibrium of 
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the asperity. This gives us c = aZmmlK. Then, du/dz = {azmm/ Kz). We 
may find the elastic energy of an asperity from 

Eelast = (1/2) / dz2'KRzK{ — f = 

[7:R{aZrmnf/K] / - = [aR{aZrmnf/K]ln{^). {AA) 

In order to compare this with the expression used in the scahng treatment of 
multiscale asperities used in Ref. [7] and here, let us now write this expression 
in terms of Am, given by 

Am = (aZmin/K) j — = {aZmin/K)ln{^^^). {5 A) 

Solving for a in terms of Au and substituting in the above expression for 
Eelast, we obtain 

jkKR 2 KnL'^Zmax I Am ^ 
= U^)^'' = 2/n(^^ ' ^^^^ 



^max 



where L is defined by ttL^ = vrr^^^ = 2T[Rzmax- For the v}'^ order asperity, 
we identify L with L„ and z^ax — Zmin with L'^. Since for all except for the 
extremes of complete compression and nonexistent compression, Zmax and 
Zmax — Zmin are of Comparable magnitudes and ln{zmax I Zmin) contributes 
only a factor of order unity, we conclude that our original scaling treatment 
of multiscale asperity multistability, which does not consider the actual shape 
of the asperities, will give correct orders of magnitude. 

7 Appendix B: Demonstration of the Occur- 
rence of Completely Plastic Contact at Very 
Small Length Scales 

In Appendix C of Persson's Journal of Chemical Physics paper on contact 
mechanics[9], he obtains a differential equation for the area of contact due 
to plastic deformation alone by integrating his diffusion equation: 

—d^-^^^^^^ ^^^^ 
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over cr, to obtain 



9cio ^^^(^'0 = /(C)[ ^U=o + ^^|.=..] (25). 

The integral Jq^ daP{a,() is by definition P{C) = Aei{C)/A. The first term 
on the right hand side of Eq. (2B) was shown in Ref. [9] to be equal 
to —dAnoniC) / dC, the negative of the derivative of the surface area not in 
contact with respect to C- In the limit as ay approaches infinity, the second 
term on the right hand side of Eq. (2B) vanishes. Since in this limit there 
is clearly only elastic deformation, the latter term must be associated with 
plastic deformation. Then in this limit, Eq. (2B) can be written as 

dA^ dAu>n ^ Q ^3^^ 

dC d( ^ ^ 

and hence Aei + Anon must equal a constant A in this limit, and that constant 
must be A. Then, clearly when ay is finite, we must identify the last term on 
the right hand side of Eq. (2B) with dApi/d(. Using his solution for P{a, () 
Persson finds by integrating this expression that 

-i2/n)±i-irC-^{l-exp[-al f' dC'giC)]}, (45) 



where a„ = nnao/ay, and 



A 



n=l ^ 



A 



-{2M±{-m'^){l-exp[-al f'dCgiC)]}- (55) 
^1 n Ji 
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Since as ( approaches infinity d('g{(') approaches infinity, we find that in 
the infinite ( hmit, the sum of Eqs. (4B) and (5B) becomes 

n=l,odd 

by a well known identity, which implies that none of the asperities at the 
regions of contact for ( = 00 are in elastic contact. Since there is no stress 
acting on the area Anon by definition and a = cry at all areas of plastic 
contact, we must have 

ayApi = aoA. {73) 

This result also follows if we take the C, approaches infinity limit in Eq. (4B), 
making use of the identity 

-x = - V(-l)"sm(n7rx), (85) 

for — 1 < X < 1. 

8 Appendix C: A Demonstration that Ar- 
chard's Model of Multiscale Asperities Pre- 
dicts that the Contact Area Per Asperity 
Increases with Increasing Load 

Archard[T6] presented a model for elastic distortion of multiscale asperities 
in which at the longest length scale an asperity has a hemispheric shape and 
distorts according to Hertz theorypH]. The area of contact of such an asperity 
with a fiat surface is assumed to be covered by a continuous distribution 
of smaller scale hemisherically shaped asperities. The largest length scale 
asperity described above is only in contact at the areas of contact of these 
smaller asperities. The smaller asperities are in turn also assumed to be 
in contact only at the locations of a continuous distribution of still smaller 
asperities. Archard argues that as we consider shorter and shorter length 
scales, the contact area becomes more and more nearly proportional to the 
load carried by the original longest length scale asperity. 

At his longest length scale, Archard finds that the contact area is propor- 
tional to W"^^^, where W is the load and the asperity density is independent 
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of load. At the next smaller length scale considered by him, the contact area 
is proportional to W^^^ and the number density of these asperities is propor- 
tional to W^^/'^, and hence the contact area per asperity (i.e., the quotient of 
these two) is proportional to W"^^^. At the next smaller length scale, the con- 
tact area is proportional to W'^^l'^'^ and the asperity density is proportional to 
VF^/^, leading to a contact area per asperity (the quotient of the contact area 
and asperity density) of M^^/^'^. Next order gives a contact area proportional 
to W'^^l'^'^ with an asperity distribution proportional to W'^^^'^'^ , leading to 
a contact area per asperity of W'^^^^^. As we continue this process, the as- 
perity density is always proportional to the same power of W as the contact 
area at the previous length scale. Hence, it is clear that as we go to succes- 
sively smaller length scales, the contact area per asperity is proportional to 
a smaller and smaller power of W. Thus, it is clear that as W increases, the 
smaller length scale asperities' contact area will grow, suggesting that the 
smallest length scale asperities will flatten out as the load increases. Since 
the amount that asperities on smaller and smaller length scales flatten out 
as a function of W decreases as we go to smaller and smaller length scales. 
Since this theory does not assign a precise radius to each length scale asper- 
ity, it is difficult to assess the degree to which smaller length scale asperities 
flatten out (i.e., how the contact area per asperity compares to the original 
size of an asperity. Since the contact area at each length scale appears to 
be proportional to the inverse asperity density raised to the same power as 
the power of W that the contact area is proportional to at that length scale, 
we can take that quantity as related to the original asperity size. Since the 
contact area per asperity grows as W to a smaller power, it is possible that 
the asperities do not flatten out in Archard's model, in the sense that the 
contact area remains small compared to the asperity size as W increases. 
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